Pattern formation in mixtures of ultracold atoms in optical lattices 
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Regular pattern formation is ubiquitous in nature; it occurs in biological, physical, and materials 
science systems. Here we propose a set of experiments with ultracold atoms that show how to 
examine different types of pattern formation. In particular, we show how one can see the analog of 
labyrinthine patterns (so-called quantum emulsions) in mixtures of light and heavy atoms (that tend 
to phase separate) by tuning the trap potential and we show how complex geometrically ordered 
patterns emerge (when the mixtures do not phase separate), which could be employed for low- 
temperature thermometry. The complex physical mechanisms for the pattern formation at zero 
temperature are understood within a theoretical analysis called the local density approximation. 



We choose to examine a physical realization of 
mixtures of ultracold atoms that illustrate pattern 
formation — a mixture of fully polarized light fermionic 
atoms with heavy atoms that are either fully polarized 
fermions or (nearly) hard core bosons. This choice im- 
plies that there are only two possibilities for the occupa- 
tion of each atomic species on an optical lattice site: ei- 
ther no particle or one particle of each species is allowed. 
We also assume that the heavy particles are essentially 
localized, but the theoretical system nevertheless has the 
freedom to rearrange itself amongst the different possi- 
ble configurations of the heavy particles (which is anal- 
ogous to the Ising model of magnetism, where the sys- 
tem does not allow dynamic quantum-mechanical spin 
flips but different spin configurations are sampled in a 
statistical-mechanical sense via an averaging over all pos- 
sible states). Most importantly, there is an interspecies 
interaction when both atoms lie on the same optical lat- 
tice site; this interaction is ultimately responsible for the 
pattern formation that emerges in these systems (pat- 
tern formation will occur for either attractive or repulsive 
interspecies interaction which are related via a particle- 
hole transformation in the homogeneous phase diagram; 
we focus on repulsive interactions here — note that attrac- 
tive interactions may have technological use in creating 
polar molecules on a lattice). Experimental realizations 
of such a system can be found, for example, in mixtures 
of heavy fermionic atoms like 87 Sr, 171 Yb, or 173 Yb with 
light fermionic atoms like 6 Li or 40 K. 

The light atoms can hop between nearest neighbors on 



the optical lattice, and they feel an inhomogeneous trap- 
ping potential (assumed to be harmonic here) that keeps 
them close to the center of the lattice. These atoms will 
be described by spinless fermionic creation cj and anni- 
hilation Cj operators associated with each lattice site i 
(because they are fully polarized). The heavy atoms do 
not hop, so they are described by classical occupation 
numbers Wi which equal zero or one at each lattice site. 
The heavy particles are in an inhomogeneous trapping 
potential (also assumed to be harmonic, but with an in- 
dependent trap frequency). Finally, there is an on-site 
interaction between the two atoms denoted by U. The 
Hamiltonian then becomes 

h = - j $>t Cj + 4 c i) + E( y * - ti°\ c i 

(i,3) i 
i i 

which is called the spinless Falicov-Kimball model [l| and 
has been proposed within the ultracold atom context 
by Ates and Ziegler ||. The symbol J is the nearest- 
neighbor hopping integral, and the angle brackets denote 
a summation over nearest neighbor pairs, while the sym- 
bols Vi [V i ) are the harmonic trapping potentials for 
the light (heavy) atoms, and /i [fi ) are the correspond- 
ing chemical potentials. For concreteness, we work on a 
two-dimensional square lattice, with lattice spacing a. 

The spinless Falicov-Kimball model is known to have 
a rich phase diagram in the homogeneous case (when the 
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trapping potentials are set equal to zero). On a square 
lattice, when the light and the heavy atoms each occupy 
one half of the lattice sites, the system is known to form 
a checkerboard pattern that alternates between the light 
and the heavy atoms for all values of U [3j, |j] . When the 
interaction is large and repulsive U — > +00, the system 
phase separates (all light atoms move to one side of the 
lattice, and the heavies move to the other side) whenever 
the total number of atoms is less than the number of lat- 
tice sites @. Elsewhere, the system displays a complex 
phase diagram [||, which includes stripe-like phases and 
more complex two-dimensional patterns Q. We expect 
the inhomogeneous system in a trap will display different 
phases at different spatial locations as the local chemi- 
cal potential sweeps through the phase diagram. This 
should produce a "generalized wedding-cake-like" struc- 
ture but with the different "layers" corresponding 
to different complex ordered patterns. 

We use a separate harmonic trap for each of the two 
atomic species (whose minima are both located at the ori- 
gin in space), and treat the trap frequency as a tunable 
parameter; this can be realized experimentally by using 
optical generated traps. By making one trap sharper, we 
squeeze the particles to lie in the center, and we can tune 
through a number of different stable phases. The tem- 
perature at which the ordering occurs also changes with 
the changing patterns, so it is possible that these studies 
could lead to interesting ways to perform thermometry. 
The phases that appear also depend on the ratio of the 
interparticle interaction U to the intersite hopping J of 
the light atoms. The ratio of these two can be tuned 
either via an interspecies Feshbach resonance that can 
change U, or by tuning the magnitude of J. 

In our calculations, we take a square lattice that is a 
50 x 50 cluster occupied with 625 light (mobile) and 625 
heavy (localized) atoms. The light atoms are confined in 
a trap potential of the following form 



Vi = J 



huj 

2Ja 



(2) 



where the trap frequency u> determines the curvature of 
the trap potential [site i has position {xi,yi)] and we 
will vary 1/30 < Tiuj/2J < 1/12.9. The heavy atoms 
are trapped in a similar potential with a fixed curvature 
corresponding to hu h /2J = 1/30. 

Our analysis is based on a quantum Monte Carlo 
(QMC) approach and on a local density approxima- 
tion [l(| (LDA) approach. The QMC algorithm is identi- 
cal to that in the homogeneous case [ll], LL2| except that 
the hopping matrix is supplemented with a diagonal term 
associated with the trapping potential. The LDA, at 
T = 0, can be constructed from the homogeneous grand 
canonical ground-state phase diagram [6j as a function of 
the chemical potentials. We determine the local chemi- 
cal potentials for each lattice site by subtracting the trap 
potential at that lattice site, and then we map out the 
quantum-mechanical phase diagram for each spatial coor- 
dinate within the trap. The chemical potentials are then 



adjusted to produce the correct total number of heavy 
and light atoms in the lattice. 

There are a number of ways to create an experimen- 
tal system that behaves like the spinless Falicov-Kimball 
model. First, one can slowly adjust the depth of the opti- 
cal lattice for the heavy atoms, to allow different config- 
urations of the heavy atoms to be sampled as the lattice 
turns on, and then freeze in one of those configurations. 
Second, one can work with heavy atoms that are not 
fully localized, but do hop slow enough that it is safe to 
ignore the quantum-mechanical energy associated with 
their motion through the lattice. Third, one can view 
each experimental configuration of the localized atoms 
as a frozen snapshot chosen according to the equilibrium 
probability distribution of the Falicov-Kimball model, so 
averaging over different experimental runs will produce 
the (equilibrium) results of the model. An imaging exper- 
iment can then measure the configuration of the atoms 
at a specific moment in time. We simulate this type 
of experimental measurement by taking representative 
snapshots of the current QMC configuration (after ther- 
malization), to see what a "typical" density profile is for 
each atomic species. 

We first show results of these snapshots for the case 
where U = 5 J (see Fig. [IJ, which is in the moderately 
strong interaction regime (the noninteracting bandwidth 
on a square lattice with nearest-neighbor hopping is 8J), 
where the system has a significant tendency to phase 
separate. We see the character of the phase separation 
change from having the heavy particles at the center, to 
having them at the periphery of the trap, as is expected 
for two species that phase separate and lie in different 
curvature traps. What is more interesting are the pat- 
terns that emerge in the QMC snapshots in the crossover 
regime, where the transition from the heavies being at 
the center to being in the periphery occurs through a 
complicated intermediate state that displays an obvious 
analog to labyrinthine pattern formation in immiscible 
ferrofluids [l3[ . Indeed, our LDA analysis shows that the 
ordered phases that create the labyrinthine patterns are 
the well-known axial stripe phases that appear over a 
wide range of the homogeneous phase diagram. Stripes 
along one axis merge with stripes oriented along a per- 
pendicular axis to create the labyrinthine effect. This 
behavior is also similar to the quantum emulsions seen in 
mixtures of different mass bosons in one dimension [14| . 

Next we examine the weakly correlated regime with 
U = J in Fig. [21 In this case, we tend to see more 
different ordered phases in addition to the phase sepa- 
ration. The system also displays a strong sensitivity of 
the pattern shape and physical extent on the tempera- 
ture (particularly the heavy atoms), which suggests that 
this could be used for thermometry purposes if the size 
of the atomic cloud could be measured with a precision 
on the order of a lattice site. One can also try to fit 
the tail of the light atom distribution (where there are 
no heavy atoms) to determine the temperature, as was 
recently done in experiments on polarized mixtures of 
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spin-one- half fermions [151 ]. 

In the panels of Fig. [51 we show representative QMC 
snapshots again, this time for only two cases — hcu/2J = 
1/17 and Tllo/2J = 1/12.9. For each case, two tempera- 
tures are plotted. One is relatively high as compared to 
the overall ordering temperature, while one is a very low 
temperature (essentially T = 0; our results are virtually 
identical for k B T = 0.0001J and k B T = 0.001J). One of 
the first things to notice is that the temperature can be 
roughly found just by examining the physical size of the 
heavy atomic cloud. With the same trap potential, the 
extent of the peripheral region of the heavy atom cloud 
depends strongly on the temperature, but needs to be 
calibrated for experimental use. The shape of the den- 
sity profile at the edges of the trap also display strong 
sensitivity to T, and might be useful for thermometry as 
well — note the ordered phases that appear at the edges 
of both low temperature panels (b) and (d) are higher- 
period ordered phases whose appearance is sensitive to 
the T. 

Finally, we compare the QMC results to the LDA at 
T = in Figs. Eland [U We show plots of the LDA for the 
density as a function of radial distance from the center 
of the trap for the QMC snapshots plotted above. We 
compare the LDA directly to the QMC average. Note 
how similar the LDA results are with those of the QMC 
simulations, especially in the cases when we do not have 
labyrinthine patterns (the difference is probably a tem- 
perature effect); the deviations for small radii arise from 
vacancies in the given snapshots. Also note that when 
there are ordered phases (and different traps), the LDA 
need not show a monotonic dependence of the density on 
the distance from the center of the trap (although the 
light atom density is always piecewise monotonic when 
the system is in one phase). 

The most direct way to observe the pattern formation 
is a tomography experiment with a resolution on the or- 
der of a single lattice spacing, which would provide an 
in situ means to directly observe the patterns. Unfor- 
tunately, such resolution is not yet available. Another 
approach is to examine Bragg scattering peaks. In par- 
ticular, when we have labyrinthine patterns forming, the 



structure factor will be peaked around (tt,0) and (0, 7r), 
while the checkerboard phase will show peaks at (tt,tt). 
The large growth of these peaks at low T determine the 
relative area of the ordered phases and provide another 
means for thermometry of the system. Finally one can 
examine the density distributions in momentum space in 
time-of-flight imaging that correlate the shape of n(k) 
with different spatial extents of the light atom clouds 
or employ noise correlation spectroscopy in time-of-flight 
imaging to examine the density-density correlation func- 
tions, whose momentum dependence will show peaks at 
wavevectors associated with particular ordering patterns. 
The momentum distribution itself (not shown here) does 
not readily show which phases are present, but it can de- 
termine the overall spatial extent of the clouds of atoms. 

In conclusion, we have shown that interesting density- 
wave patterns appear in mixtures of spin polarized 
fermionic atoms with different masses. The ordering phe- 
nomena ranges from geometrical phases to labyrinthine 
patterns reminiscent of viscous fingering in immiscible 
ferrofmids. The ordering phenomena is most striking 
when the heavy mass atom is localized on the lattice, 
but we expect these patterns to survive when the kinetic 
energy for the heavy atoms remains small (due to a Fermi 
pressure effect [111). Note that one needs to achieve low 
temperatures to see these patterns. While this is not 
yet feasible in current experiments, new developments in 
cooling technology [l7| make it seem likely that this phe- 
nomena will be able to be observed soon. 
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FIG. 1: (Color online) Representative QMC snapshots of real-space configurations of heavy (top) and light (bottom) atoms, 
respectively. All simulations have been performed for U = 5J at temperature feT = 0.01J. Each image pair as we move 
from left to right corresponds to a different huj/2J value: (a) 1/30; (b) 1/20; (c) 1/18.5; (d) 1/17; and (e) 1/12.9. The heavy 
atoms in a snapshot either occupy a site (light gray or red) or do not (black), while the light atom densities are plotted with a 
false-color plot. The heavies migrate from the center to the periphery as the light trap tightens (left to right). 




FIG. 2: (Color online) Representative QMC snapshots of real-space configurations of heavy (top) and light (bottom) atoms, 
respectively. These simulations were performed for U = J. The parameters are as follows: (a) k B T — 0.05J, hu>/2J — 1/17; 
(b) k B T = 0.0001J, huj/2J = 1/17; (c) k B T = 0.05J, hu/2J = 1/12.9; (d) k B T = 0.0001J, hu/2J = 1/12.9. 
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FIG. 3: (Color online) Averaged heavy and light radial densities per lattice site for U — 5J and ksT = in the LDA (dotted 
line) compared to the QMC snapshot average at ksT — 0.01J (solid line). The five cases are identical to those shown in Fig.[T] 
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FIG. 4: (Color online) Averaged heavy and light radial densities per lattice site for U — J and ksT — in the LDA (dotted 
line) compared to the QMC snapshot average at feT = 0.000IJ (solid line). The two cases are (a) Ttu)/2J = 1/17 and (b) 
hu/2J — 1/12.9 and correspond to the low T images in Fig. [2] 



